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I. INTRODUCTION 
Hydrodynamic equations 


the values of second-order terms and beyond become 


negligible and can be disregarded. 


There is a lack of prior research systematically investigating 


are commonly expressed 
through a truncated Taylor series, where terms of order 2 
and higher are omitted, reducing the series to include only 
order 1 differentials. This simplification is based on the 
assumption that, within small intervals such as, ôt, ôx, 62, 


Abstract— The Weighted Taylor series is an adaptation of the 
conventional Taylor series truncated to the first order, wherein high- 
order differential terms are replaced by introducing weighting 
coefficients to the initial terms. This study presents a novel approach 
for computing these weighting coefficients specifically designed for 
water wave modeling. Subsequently, the derived weighted Taylor 
series is employed to formulate both weighted continuity and the 
weighted Laplace equation. The weighted Laplace equation facilitates 
the formulation of the velocity potential equation, leading to the 
development of wave 
important phenomena such as shoaling, breaking, and refraction- 
diffraction. Additionally, the formulation of the weighted Euler 
momentum conservation equation is introduced to determine the wave 


transformation equations encompassing 


number in deep water. By scrutinizing the outcomes of dispersion 
equations, as well as analyzing shoaling-breaking and refraction- 
diffraction scenarios, optimal values for the weighting coefficients are 
identified.. 


weighting coefficient in the truncated Taylor series can alter 
wave characteristics, shortening wavelength and reducing 
water particle velocity. These findings suggest that utilizing 
a weighted Taylor series with optimized coefficients can 
enhance water wave models. 


In this study, the Taylor series was truncated using the 
central-difference method. This method only extracts 
contributions from even-order differential terms, and is 
subsequently corrected by the remaining contributions from 
odd-order differential terms. 


the Taylor series truncation method. Hutahaean (2021) 


introduced the weighted Taylor series by examining 
intervals where the values of 1st order terms significantly 
surpass those of 2nd order terms. Building upon this, 
Hutahaean (2022) applied the Forward Difference scheme 
to formulate the weighted Taylor series, assigning a 
weighting coefficient solely to the time derivative term in 
the function f(x, z, t). Both studies demonstrated that the 
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The truncated Taylor series was then utilized to formulate 
the foundational equation of hydrodynamics, specifically 
the continuity equation. Each term in this equation received 
its respective weighting coefficient, leading to the 
designation of a weighted continuity equation. 


Employing the weighted continuity equation, the 
formulation of the weighted Laplace Equation followed. 
Subsequently, the weighted Laplace equation was solved 
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for sloping bottoms to derive conservation equations 
governing variations in wave constants as waves transition 
from deeper to shallower waters. 


By utilizing these conservation equations, the study derived 
equations describing changes in wave constants, including 
wave amplitude and wavelength in shallower waters. These 
equations incorporate weighting coefficients, with the 
corresponding values of the weighting coefficients studied 
based on the results of dispersion equation, shoaling 
breaking model, and refraction-diffraction model. 


II. WEIGHTED TAYLOR SERIES 
Taylor series for a function with two variables, f = 
f(x,t), (Arden, Bruce W. and Astill Kenneth N. ,1970) is 


f(x + ôx, t + dt) = f(x, t) + +6¢ + 


a ôt? a? a Sx? a? 
gu Ee See a 
ax 2! at? dtdx 2! dx? 


(1) 

In this paper, x represents the horizontal axis, z the vertical 
axis, and t time. In the formulation of hydrodynamic 
equations, including continuity equation and Euler's 
momentum conservation equation, this Taylor series is 
truncated to a first-order differential series, assuming that 
within the interval dt, ôx are very small, higher-order 
differential terms are very small and can be neglected, i.e., 
it becomes:. 


df df 
f(x + 6x,z + 6z,t + ôt) = f(x,z,t) + état Ox a 
However, as the interval decreases, not only do the values 
of higher-order differential terms decrease, but the values of 
first-order terms also decrease. As a result, the relative 
values of higher-order terms to the first-order terms remain 
significant. 
Equation (1) can be written as follows, 


f(x + 6x,t + 6t) = f(x,t) + 


teen ) oF 
( 2a °” at 
Howe = )s a 
( 2 ay) ae 


There are contributions from higher-order terms to the first- 
order terms. Therefore, these higher-order terms cannot be 
simply neglected. Truncating the series to the first order 
should be accompanied by providing coefficients that 
represent the higher-order terms. 

f(x + 6x,z + 6z,t + ôt) = f(x,z,t) 


df df 
+ y.6t = + Vx 5X ax 
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This series is a weighted Taylor series, with weighting 
coefficients y, and yy. With this equation, the role of higher- 
order terms is not completely eliminated; it is represented 
by the weighting coefficients. The characteristics of the 
function present in the higher-order terms are still reflected 
in the weighting coefficients. 


The aim of this research is to develop a method for 
calculating weighting coefficients by extracting 
contributions from higher-order terms and adding them to 


the first-order terms. 


The method of absorbing contributions from high-order 
terms or the formation of weighting coefficients generally 
consists of two parts: 


a. Absorbing contributions from odd-order 
differential terms. 
b. Absorbing contributions from — even-order 


differential terms. 


In this research, weighting coefficients will be formulated 
for the water wave modeling equation. The solution to the 
velocity potential equation of the Laplace equation (Dean 
(1991)) is: 

o(x,z,t) = G cos kx cosh k(h + z) sin øt 


: 2T . 
k is wave number, k = 7 ois the angular frequency, o = 


2 f : ; ee 
E T is wave period and h is water depth. Considering that 


sin function has similar characteristics to the cos function, 
a method for calculating weighting coefficients will be 
developed using a functional form., 

f(x,z,t) = cos kx cosh k(h + z) cos øt 


MI. ABSORPTION OF CONTRIBUTIONS 
FROM THIRD-ORDER DIFFERENTIAL TERMS 
In this research, a calculation method is developed 

wherein the Taylor series is employed up to the third order 
only. Limitations are introduced by utilizing an interval size 
where terms of the fourth order and higher become 
negligible. Subsequently, the third-order terms are removed 
from the series, and their values are added to the first-order 
terms. The magnitude of the contribution of third-order 
terms to the first-order terms is expressed with a 
contribution coefficient. 
The contribution coefficient for the time differential terms 


is given by, 
ôt? arf 
_ 6 dt? 
Ht = spat 
dt 


As a function of time t is f(t) = cos øt, the differential 
substitution of this function, performed under the condition 
cos ot = sin ot 
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2 
=—-—oO 
Ht 6 
Substituting ôt = ¢,T and o = 2 
T 
An? 
hmn 


€ is the time interval coefficient t. Similarly, for the 
function f(x) = cos kx, the following is obtained 


4r? 3 
Ux = -5 & 
€, is the coefficient of the interval x, where 6x = €,L, and 
Ar? 5 
Uz = p? 


€, is the coefficient of the interval z, where 6z = €,L. 


It is important to note that the computation of contribution 
coefficients necessitates values for the interval coefficients 
Et, Ey and €, which will be discussed in the subsequent 
section. 


3.1. Calculation of Interval Coefficients 
In this section, the method of calculating the values of 
interval coefficients €+, €y and €z, is discussed. The 


calculation is performed using the optimization equation. 
S2+s3 


S1 


S1, S2 dan S3, are, in order, the first, second, and third terms 


<e 


in the Taylor series. Term 1 is a term with a first-order 

differential, term 2 is a term with a second-order 

differential, and term 3 is a term with a third-order 

differential. 

€ is a small number whose value is determined. The larger 

the value of £, the larger the values of higher-order terms 

that will be extracted. Furthermore, € is referred to as the 

optimization coefficient.. 

The optimization equation is an equation with variables 

representing the size of the interval, i.e., ôt, ôx, ôz. This 

equation for the interval size is converted into an equation 

for the interval coefficients, i.e., &&, €, and €z, by 

substituting ôt = ¢,T, dx = €L dan 6z = €,L. T is the 

wave period, and L is the wavelength. 

a. Calculation of Time Interval Coefficient at -t £, in the 
function f(t) = cosot. 

By utilizing a series up to order 3, the optimization equation 


can be expressed as, 
strap dt3df 
2 dt? 6 dt3 


=E 
af 
ota. 


Substituting the differentials of f(t) when cos ot = sin ot, 
a second-degree polynomial for St is obtained. 

ôt ôt? , 

=o —— 0t =E 


2 6 
Substituting ôt = ¢,T , and o = = where T is wave period, 
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2m? , 
z éi —Té& te=0 


Of the two €+, the least value is used. 


In the equation for €; there are neither variable o or T. 
Hence, € is not dependent to the wave period T and it 
applies to any wave period T. 


b. Interval Coefficient -x £x on function f(x,t) = 
cos ot cos kx. 

For this function, there will be two variables in the 
optimization equation €, and €. By using € from the 
calculation results for the function f (t) = cos ø t, hence £x 
is the only unknown variable. 

The optimization equation for the function f (x, t) is worked 
out using the Taylor series up to the third term, resulting in 
a third-degree polynomial. This polynomial can be solved 
using the Newton-Rhapson iteration method, which 
requires an initial iteration value. To obtain the initial 
iteration value, the optimization equation is first solved 
using only the second-order differential terms, forming a 
quadratic polynomial.. 


23 >e 
S1 
Where 
f df 
S1 = état xa 
= ot? a? f Zf = bx? df 


saa t 6tdx ==> 
2 dt? dtdx 2 dx? 
Substituting s, and s, into the optimization equation, the 


S} = 


computation proceeds under the conditions cos ot = sinot 
and cos kx = sin kx, resulting in the cancellation of terms 
in the numerator and denominator. Subsequently, the 
substitution ôt = ¢,T, dx = €,L,0 = = and k= = with Et 
as a known variable, a quadratic equation in €, is derived 
from calculations involving the function f(t) = cos ot 
Cot Cex + C E = 

Co = 2TE,E — 2m" Ef 

Ci = 2me + 4n? E, 

C3 = —2n” 

The quadratic equation yields two roots, of which the larger 
one is selected. Notably, the €, \ equation no longer contains 
the wave period T or wavelength L. 


Proceeding to the optimization equation (2), we obtain a 
cubic polynomial. 


ôt’ df ôt? d3f ôx? af 
S3 = + ôx 

6 dt? 2 dt?dx 2 dtdx? 

ôx? def 

6 dx? 
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Where, 

Co + Cex + C2 EŻ +063 62 = 0 
813 

Co = 2mepe — 2T? Ef + re 

2 81° 2 

C1 = ZTE + 4N Et t- ét 


5 8r? 
C = —2rm^ + t 
8r? 


The solution to this equation can be determined using the 
Newton-Rhapson iteration method, initiating the process 
with the root derived from the quadratic equation. 


c. Interval coefficient -z €, for the function f(x, z, t) = 
cos øt cos kx cosh k(h + z). 
The equation for calculating ¢, is formulated using a 
method similar to the equations used to calculate ¢, and €,, 
namely by utilizing optimization equations. By applying 
conditions such as ccosot = sin øt, cos kx = sin kx dan 
cosh k(h +z) = sinh k(h + z), specifically when h is 
large and z is very small, the elements in the numerator and 
denominator cancel each other out. Substituting dt = ¢,L, 


21 2T : 
ôx = E&L, 6Z = €,L, o = T and k = 7 we obtain an 
equation with three variables: €,, €, and €,, with €, and Ey 


obtained from previous calculations, one variable, and ¢, 
remains. 


The quadratic equation for €, is, 

Co2 + C12& + C2267 = 0 

Coz = (& + €,)2ue — 2072? + 4n?’ EEx — 20782 
C12 = —2 me — 4m E, — 4n?’ E, 


C22 = 21? 


This quadratic equation yields two values for €,, with the 
larger value chosen for subsequent analysis. These solutions 
form the input for a third-degree equation. 


813 5 8r? , 813 3 813 7 

Co3 = A E” + 5 Ex Ep + 5 EtEx + eS Ex 
Co = Co, + Co,2 

8r? 7 , 8r? . 
C13 = = êt + 81° Ep Ey ara & 
c1 = 2+ C13 

813 813 
C23 = aS Et — 2 Ex 


C2 = C22 + C23 
873 
C3 = — 


6 


Co,2» €1,2 and C2 2 are derived from the quadratic equations.. 
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Table (1) below showcases examples of interval coefficient 
calculation results for various optimization coefficient 
value € . 

Table 1. The values of interval coefficients 


E Et Ex E, 
0.005 0.0016 | 0.00478 | 0.01452 
0.006 | 0.00192 | 0.00574 | 0.01747 
0.007 | 0.00224 0.0067 | 0.02044 
0.008 | 0.00256 | 0.00766 | 0.02342 
0.009 | 0.00288 | 0.00862 | 0.02641 
0.01 | 0.00321 | 0.00959 | 0.02942 


As seen in Table 1, as the value of € increases, the interval 
coefficient also increases, indicating a greater contribution 
from higher-order terms. With these interval coefficients, 
the contribution coefficients of the higher-order terms are 
calculated. 

Table 2: The values of the contribution coefficients. 


E Me Hx Mz 
0.005 | -0.000017 | -0.000151 | 0.001388 
0.006 | -0.000024 | -0.000217 | 0.002009 
0.007 | -0.000033 | -0.000296 | 0.002750 
0.008 | -0.000043 | -0.000386 | 0.003611 
0.009 | -0.000055 | -0.000489 | 0.004596 
0.010 | -0.000068 | -0.000604 | 0.005705 
IV. CONTRIBUTION FROM EVEN 

DIFFERENTIAL TERMS 


In this stage, the Taylor series is calculated with a very 
small interval size, so the fourth-order terms and higher are 
very small, while the third-order term has been absorbed by 
adding its value to the first-order term. Thus, the remaining 
Taylor series of the second order takes the form, 


f@+6x%,t+6t)=f@zon+ + uote + 


1+ Uy). Lor aie 1 
rhe ag oe ap tag 
bx? d'f 
2 ae 
The simplification of the subsequent expressions is as 
follows. 
a, = 1+ 
a, =1+ ux 


Truncated Taylor series becomes 


d 
f(x + 6x,t + ôt) = f(x,z,t) + asta 
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tages ee ates ee 
2 dt? dtdx 2 dx? 
This Dn ta involves the absorption of contributions 
from even higher-order differential terms through the 
application of the Central Different Method operation. 
4.1. Function f (x, t) 


Truncated Taylor series for (x + 6x, t + ôt) 


d 
f(x + 6x,t + 8t) = f(x,z,t) + at L4 
af ôt? df Zf 8x af 


Be oh a 
oe ded Ty T 72 Ie XJtax 2 dx? 


Truncated Taylor series for (x — 8x, t — ôt) 


d 
f(x — 6x,t — ôt) = f (x,z,t)— ast L- 
af ôt?a? d? ôx? d? 
ay ôx Hoa ôx f r 
dx 2 dt? dtdx 2 dx? 


These two equations are subtracted from each other, 
f(x + 6x,t + 8t) — f(x — x,t — ôt) 
df df 


= 2a,6t — + 2a,6x — 
a ar t a,0x = 


During this operation, the even higher-order differential 


; ri af i a 
terms automatically vanish. 3 will be absorbed at s, and 


the last equation is divided by 26x, 
f(x + 6x,t + dt) — f(x — 6x,t — ôt) 
26x 


_ ot df : df 

= t xat ax 
For very small values of ôt and ôx, this equation becomes 
a total spatial differential, 
D f ôt d T; df 
dx “tsx dt tax dx 
The Taylor series, which has eliminated both even and odd- 
order differential terms, is, 


d 
f(x + 6x,t + ôt) = f(x,z,t) + mot SL 
af 
y- 
+a: Xay 
Substituting the total spatial differential into -a 
f(x + 6x,t + 6t) = f(x,t)+ a,0t r af zt Qx 
ôt df df 
(a Sx ae + a,,0x =) 
df 
fx +6x,t+ 6t)=f@w%o+d+ Ox) Opt = 
df 
265 sie 
+a% Xas 


This equation is written as, 


f(x + 6x,t + ôt) = f(x,t) + yt2ôt oxo df 


www.ijaers.com 


International Journal of Advanced Engineering Research and Science, 10(12)-2023 


This equation represents a weighted Taylor series for the 
functionf = f (x,t) with weighting coefficients, 

Ye2 = 1+ ay) a; 

Vx = az 

The use of the index 2 in y;, indicates that the weighting 
coefficient is for a function with two variables, f = f (x,t). 
The cross-contribution, i.e., the contribution of spatial 
derivative to the time derivative, is present in @,,@;. In cases 
where contribution coefficients are negligible, where yu, and 
Lx are both zero, resulting in @, = a, = 1, we obtainy,, = 
2 and y,=1, these values serve as the fundamental 
weighting coefficients for the function f (x, t). 


4.2. Function with three variables f = f (x, z,t). 
To obtain the contribution of even-order differential terms 
for a function with three variable f = f(x, z,t), the total 


: aoaea D ; : ; 
spatial derivative 21 from the function with two variables 
f =f(,t) and 21 on the function f = f (z,t) where the 
formulation of 21 is performed in the same method as the 


i D aa 
formulation of or obtaining 


Df _ sot af af 
Ke ~ wt Brae taza 
ea of and 2% is substituted into Ê = L and = = Lon the Taylor Series 


f=f te z,t) which has seen its higher-order terms, 
we obtain 


d 
f(x + 6x,z + 6z,t + ôt) = f(x,z,t) + Te 


d 
+y,ôx Fy, 82 z 


Yt = (1+ a, + @,)a, 
Vx = az 
Yn až 


It is observed that the value of yy on f(x, z, t) is the same 
as the value of yy for the function f(x, t). In cases where 
contribution coefficients are neglected, where py = Uy, = 
Uz = 0 resulting in a, = a, = a, = 1, we obtain y,3 = 3, 
while yy = y, = 1. These values serve as the fundamental 
weighting coefficients for the function f (x, z, t). 

The values of weighting coefficients for various values of € 
are presented in Table (3). 


Table 3:. The values of weighting coefficients. 


E Yt,2 Yt3 Yx Yz 
0.005 | 1.99990 | 3.00237 0.9994 | 1.00556 
0.006 | 1.99986 | 3.00344 | 0.99913 | 1.00805 
0.007 | 1.99980 | 3.00471 | 0.99882 | 1.01103 
0.008 | 1.99974 | 3.00619 | 0.99846 | 1.01450 
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0.009 | 1.99967 
0.01 | 1.99960 


3.00789 
3.00980 


0.99804 
0.99758 


1.01847 
1.02295 


It is evident that the values of the weighting functions are 
not significantly different from the fundamental values, 
namely ¥,2 = 2.0, y¥z,3 = 3.0, yy =1.0 and y, = 1.0. 
However, to satisfy a balance equation, it is advisable to use 
accurate values for the weighting coefficients. 

As € increases, the values of y;2 and yy, decrease, while 
Yt and yz increase. Considering that higher-order terms 
contain information about the function's characteristics, it is 
preferable to choose € such that the coefficients. To obtain 
the optimal £ coefficient, precise calibration of the model 
results formulated using the weighted Taylor series is 
necessary. In the modeling of water wave transformation, 
this calibration can be applied to evaluate the breaker height 
or breaker depth generated by the model. 


Vs THE APPLICATION ON WATER WAVE 
MODELING 
The formulation of various equations in this paper is not 
provided in detail; only the final results are presented. The 
research focus is on the weighting coefficients. 
5.1. Weighted Laplace Equation and its Solution 


a. Weighted Laplace Equation 


The continuity equation formulated using the weighted 
Taylor series takes the form of a weighted continuity 
equation, 

du dw 
u(x,z,t) is the horizontal water particle velocity, and 
w(x,z,t) is the vertical water particle velocity. The 
weighted continuity equation cannot be expressed as 


nau | aw _ 
yzdx dz 
Nor 

du nwo 
dx y, dz 


since it yields different solutions from (3). Therefore, it can 


be said that yẹ s directly related to z and y, is directly 


a 
related to —. 
dz 


a ? i f a 
Substituting the velocity potential properties, u = — ze and 
a : : baat ; 
w=- 2, into the weighted continuity equation, the 


weighted Laplace equation is obtained, 


2 2 
ae, Foo 


Yx dx? Z ąz? 
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This equation is referred to as the weighted Laplace 
equation. 


b. Solution of the Weighted Laplace Equation 


The solution of (4), obtained using the separation of variable 
method and working with the kinematic bottom boundary 
condition for a sloping bottom, results in the velocity 
potential equation, 
p(x,z,t) = GB (z) cos kyx sin ot + 

GP (z) sink,x sin ot 


In this equation, there are two wave numbers: the wave 
number in the horizontal axis\ 


k . . : 2 
k, = re the wavelength in the horizontal axis L, = = and 
Yx x 
; P k 
the wave number in the vertical axis k, = T the 
zZz 


: : : 2 
wavelength in the vertical axis L; = = 
Zz 


ae kz(z+h) + e7kz (z+h) 


2 


kz(z+h) _ eKz(zth) 


2 


Biz) = 


ae 


B, (2) = 
i 1 dh 

Viz V¥x %* 

1 i dh 
VV JY %* 


dh. ane ; : 
aS the bottom slope, which is negative for waves moving 


a 


from deep water to shallow water.. 
On a = 1, B(z) =coshk,(h +z); 6,(z) = sinhk,(h + 
z) 


c. Equation for G 


The equation for G is obtained by integrating the weighted 
Kinematic Free Surface Boundary Condition with respect to 
time (Hutahaean (2023b), 


G= Yt 24 
~ 1 kA 
2k (== - $F) Bon) 
We 2) 
On -On 
ae" +e 
pOr) = ~z 
@ is the deep water coefficient, where man = 1. In this 


research, 0 = 1.95. In the analysis of shoaling-breaking, 
this coefficient plays a role in determining the breaker 
depth. A larger 0 corresponds to a deeper breaker depth, 
while a smaller 0 results in a shallower breaker depth. 
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d. Wave Amplitude Function 
The equation for G can be expressed as the wave amplitude 


function, 
F 2Gk Bco ( 1 -) 
= TU — 
OYt2 V¥z 2 


d. Water Surface Elevation Equation 

n(x, t) = A cos kx cos ot 

e. Wave Number Conservation Equation 

dk (n+) 
dx 

f. Energy Conservation Equation 


=0 


5.2. Dispersion Equation 


By using the weighted Taylor series, and by using the 
conditions in the formulation of the continuity equation 
where in a control volume the horizontal velocity only 
changes on the horizontal axis and the vertical velocity only 
changes on the vertical axis, and by working on Newton's 
principle of conservation of momentum we obtain two 
equations momentum, namely the momentum equation in 
the horizontal direction and the momentum equation in the 
vertical direction, respectively, are, 


du du 1 dp 
Yes ae + Yat a ~ ~ pdx 

dw dw 1dp 
Yaa, t YzaW aS = “paz 4 


The vertical momentum equation is expressed as the 
f ; . a 
integration of the pressure equation = over the water depth. 


By employing surface dynamic boundary conditions p, = 
0, the pressure equation p is obtained. Subsequently, 
differentiating the pressure equation with respect to the 
horizontal axis and substituting it into the horizontal 
momentum equation, we apply this at z = 7 to yield the 
surface momentum equation. 


du, Yx du,u dn 
my Px Senn 
Yes ae 2 ax I ax 
The dispersion equation is obtained by utilizing the surface 
momentum equation while disregarding the convective 


acceleration term, 


Substituting the velocity potential by employing its inherent 
property and substituting the water surface equation and the 
wave amplitude function, the dispersion equation is 
derived,, 
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A 
ae = = k+ Yt2Vt30° =0 


This dispersion equation is solely for the calculation of 
wave numbers in deep water. To determine wave numbers 
in shallow water, the shoaling-breaking model is utilized.. 


5.3. Shoaling Breaking Equations 


By working on the conservation equations, the G dan wave 
amplitude function (Hutahaean (...)), he shoaling-breaking 
equations are derived. For waves transitioning from x to 
x + 6x, 


dk 4k dh 

ax (4h +3A)dx 
dk 

kx4+5x = ky + Ôx 

dA G dk/1 kA 

ar (a T) 
dA 

Arra = A, T Ox a 


1 
Geiss — en Gy—5(n ky46x—In ky) 


5.4. Refraction-Diffraction Equations 


The shoaling-breaking equations can be transformed into 
refraction-diffraction equations using the procedure 
outlined by Hutahaean (2023a). 


VI. OUTCOMES OF THE MODEL 

a. The Results of Dispersion Equations 

In the following section, the computed wavelengths for 
waves with a period of 8 seconds are presented. The wave 
amplitude A varies, with the wave height H = 2A. An 
optimization coefficient of € = 0.01 where y;,. = 1.99960 
»Yt3 = 3.00980, yx = 0.99758, yz = 1.02295.. The 
calculated results are summarized in Table (4) as follows. 

Table 4: Wavelength from the modelling outcome 


H L Ly Lz H 
(m) | m (m) (m) Ly 
2 12.132 12.117 12.264 0.165 
2.1 11.784 11.77 11.912 0.178 
2.2 11.398 11.384 11.523 0.193 
2.3 10.959 10.946 11.078 0.21 
2.4 10.434 10.422 10.548 0.23 


2.5 9.739 9.727 9.845 0.257 


The calculation results for the wavelengths reveal the 
presence of two wavelengths: the horizontal wavelength Ly 
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and the vertical wavelength L,, both exhibiting small 
differences. 
Furthermore, it is observed that within one wave period, the 
wavelength decreases as the wave height increases, 
exhibiting different wave steepness. The critical wave 
steepness according to Michell (1893) is 0.142. The critical 
wave steepness according to Toffoli, A., Babanin, A., 
Onaroto, M., and Wased, T. (2010) is 0.170, potentially 
reaching 0.200. The model results closely align with the 
critical wave steepness from Toffoli et al., specifically at 
wave heights of 2.1-2.3 m. 
In conclusion, the critical wave steepness from Michell 
(1893) and Toffoli et al. (2010) indicates the maximum 
wave height for a given wave period under undisturbed 
conditions, specifically in deep water. Referring to Toffoli 
et al.'s criteria, the maximum wave height for waves with a 
period of 8.0 sec. is expected to be 2.30 m. However, in this 
study, the maximum wave height for a wave period of 8.0 
sec. is found to be 2.50 m, with a wave steepness of 0.257. 
The relationship between wave period and wave height, 
according to Wiegel (1949, 1964), is expressed as, 
gT? 

15.6? 
For waves with a period of 8.0 sec., the calculated value is 
H = 2.58 m.. In this case, the model results closely match 
the Wiegel (1949, 1964) equation with £ = 0.01. 


b. The outcomes of shoaling breaking modelx 
The shoaling-breaking model is applied to waves with an 


amplitude A = 1.20 m in coastal waters with a bottom slope 


dh ae ; 
Ta —0.005, considering various values of € and a deep 
x 


water coefficient 9 = 1.95.. The wave period is calculated 
using the Wiegel (1949, 1964) equation: 


r=156 | (sec) 
= 3 —= sec 
g 


Given a wave amplitude A = 1.20 m, H = 2.40 m, the T = 
7.716 sec is obtained. 

According to Komar, Paul D., and Gaughan, Michael K. 
(1968), the breaker height is 

Hp = 0.39 g /5(THg) /5 m. 

Hy is the deep water wave height. For waves with Hp = 
2.40 m, T = 7.716 sec the Komar-Gaughan equation yields 
H, = 2.809 m . In Table (6), on € = 0.01, H, = 2.812 m. 

Table 6:Breaking conditions for various values € 


e H, | h, | Leo | m | m, 
(m) (m) m) | Lx» hp 
0.005 | 2.883 | 3.619 | 4.539 | 0.635 | 0.797 
0.006 | 2.876 | 3.615 | 4.532 | 0.635 | 0.796 
0.007 | 2.866 | 3.607 | 4.521 | 0.634 | 0.794 
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0.008 | 2.854 3.6 | 4.509 | 0.633 | 0.793 
0.009 | 2.837 | 3.588 | 4.491 | 0.632 | 0.791 
0.01 | 2.812 | 3.568 4.46 | 0.631 | 0.788 


The results of the breaking conditions for various values of 
€ are presented in Table (6). From £ = 0.005 to € = 0.01, 
the breaking conditions vary only in the third decimal place, 
with ôe changing by 0.001. However, the difference 
between the breaking conditions with € = 0.005 and ¢ = 
0.01 is quite significant. 

The breaker depth index from the model is = = 0.788, 
which closely aligns with Mc. Cowan's (1894) criterion of 
= = 0.78. However, the most influential factors in 
determining the breaker depth and breaker depth index are 
the deep water coefficient 0, the larger 0, the smaller the 
breaker depth index. The contribution coefficient only 
corrects the third decimal place. In Table (6), the breaker 


depth index is obtained using 0 = 1.95. 


The smallest breaker steepness 2 is 0.631, achieved at € = 
ab 


0.01. This breaker steepness is significantly larger than the 
critical wave steepness from both Michell (1893) and 
Toffoli et al. (2010). This is because the breaker steepness 
occurs in regions of high wave energy concentration. 


Breaking occurs when (= — =£) = 0, 


V¥z 2 
kA 1 
2 Vy 
Considering ky, = = or k=ky,./y, and ky = = the 
following is obtained, 
Hy 2 


Lyp T [YN Yz 
Hence, the breaking characteristics are determined by the 
coefficients of the truncated Taylor series. The smaller yy 


: oe H 
and y, then the higher the critical wave steepness os 
x,b 


making breaking more challenging. Conversely, the larger 


zi H ; 
Yx and y, the smaller the critical wave steepness a making 
‘x,D 


waves more prone to breaking. 


c. Results of Refraction-Diffraction Model 

The equations derived from the shoaling-breaking 
equations can be transformed into refraction-diffraction 
equations using the method proposed by Hutahaean (2003). 
The refraction-diffraction model is executed with two € 
values: 0.005 and 0.01, considering the bathymetry of a 
submerged island (Fig (1)). Contour plots of 2D wave 
height are presented in Fig (2), with (a) using £ = 0.005 
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and (b) using £ = 0.01. Three-dimensional contour plots 


are depicted in Fig (3). 
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Fig.1: Contour of Submerged Island Bathymetry 
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Fig.2: Contour of Wave Height 


Fig.3: 3D Contour of Wave Height 
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Both refraction-diffraction model results show no 
discernible differences. From these findings, it can be 
concluded that the contribution coefficient does not play a 
significant role; the primary factors at play are the main 
values of weighting coefficients. 


In conclusion, the results of the study on weighting 
coefficients are as follows: 

The baseline values of weighting coefficients are, 

For the function f (x,t): Y2 = 2.0, Yy = 1.0 

For the function f (x, z, t): Yt3 = 3.0, yy = 1.0, y, = 1.0 
The recommended corrected values for the weighting 
coefficients are based on the optimization coefficient ¢ = 
0.01 where, Y;2 = 1.99960 ,Yt3 = 3.00980, y, = 
0.99758, y, = 1.02295. 


VIL. CONCLUSION 


As widely recognized, the precision of the Taylor series 
hinges on the inclusion of high-order terms. The greater the 
number of these terms, the more accurate the approximation 
becomes. High-order terms encapsulate essential features of 
the function under consideration in the series, making them 
indispensable. When restricting a series to only first-order 
terms, the incorporation of a set of weighting coefficients 
becomes necessary to account for the high-order terms, 
resulting in a weighted Taylor series. 


The Central Difference Method establishes the foundational 
values for these weighting coefficients. These baseline 
values undergo refinement through the inclusion of 
contribution coefficients derived from high-order terms 
with odd differentials, introducing relatively minor 
adjustments. In the context of modeling shoaling breaking, 
these contribution coefficients make subtle corrections to 
breaking characteristics, particularly at the third decimal 
place. Notably, the baseline values of weighting coefficients 
are derived without a prerequisite knowledge of the 
functional form, while determination of contribution 
coefficients relies on knowledge of the functional form. 
Thus, it is evident that baseline values hold a general 
applicability across various functional forms. 


Despite the modest impact of contribution coefficients on 
both weighting coefficients and model outcomes, their 
significance should not be understated, as they play a crucial 
role in ensuring accuracy in balance equations. In models 
demanding high numerical precision, such as time series 
models, it is advisable to utilize weighting coefficients 
corrected by contribution coefficients to enhance accuracy. 
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